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1 • intr oduc tion 



This thesis examines the use of a microcomputer to solve 
the internal processing functions of the current P-3C Update 
Omega navigation system. The most important reason for 
studying a microcomputer for the Omega navigation system is 
that during an excessive system load in the P-3C aircraft, 
the Omega signal processing is partially inhibited. 
Operating in this inhibited mode appeared unnecessary in a 
sophisticated system such as the P-3C incorporates, 
especially with the existing technology of microprocessors 
and distributed systems [Refs. 1, 2, and 3]. With a 
microcomputer dedicated to Omega navigation processing, this 
inhibited mode of operation could be eliminated ana position 
updating information could be available continuously. 

The analysis of the existing velocity and navigation 
processing routine utilizes the current equations in the 
published manual on the P-3C Update system functional 
description for Omega. The main tool for the analysis is 
the data flow graphs of the existing equations. This gives 
efficient and optimal PL/M code. 

Section II of this thesis presents fundamental 
information applicable to the worldwide Omega navigation 
system and the AN/ARN-99(V) receiver-converter. Section III 
introduces the various internal processing functions of the 
system and presents a data flow graph depicting the 
parameters which apply to each of the functions. Section IV 
presents the analysis of the velocity and navigation 
processing routine with applicable data flow graphs. 
Section V gives the conclusion of the study. 
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II. 



THEORY OF OMEGA SYSTEM OPERATION 



OMEGA is a worldwide, very low frequency (VLF) , radio 
navigation system which utilizes the radiation from eight 
transmitting stations to provide global coverage to 
aircraft, ships, land- vehicles , and submarines for accurate 
and reliable positioning. Table I lists the presently 
available transmitting stations with appropriate letter 
designators, coordinates, and direction cosine values [Ref. 
4]. The position of each transmitting station is shown in 
Figure 1 . 



A. SYSTEM PRINCIPLES 



Each of the eight transmitting stations radiates 
continuous, sinusoidal wave bursts at 10.2 kHz, 13.6 kHz, 
and 11.3 kHz. These signals are phase locked and time 
synchronized to Universal Time such that ail three signals 
start at zero value with positive slope at 0000 hours 
Greenwich Mean Time (GMT) and repeat at ten second intervals 
[Ref. 5]. Figure 2 depicts a standard ten second interval 
transmission pattern. 

In order to understand the complexity of the Omega 
navigation system and the operations performed by the 
receiver and the computer software functions, a simplified 
view of the systems operation is presented first. Assume 
some arbitrary transmitting station is radiating on one of 
the three frequencies, and, at some distance d, a receiver 
is acquiring the signal for future processing. The 
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STATION 


LOCATION 


L AT/LONG 


DIRECTION COSINE 


A 


Norway 


66 ° 25' 
13° 8 * 


12.39" N 
12.65" E 


SA 1 = +0.91 551829 
SA2 = +0.39172334 
SA3 = +0.09153733 


B 


Liberia 


6 ° 1 8 ' 
10 ° 391 


19.39" N 
44.21" W 


SB 1 = +0.98221041 
SB2 = +0.98117517 
SB3 = -0. 16588382 


C 


Ha waii 


2 io 2 4 ' 
157° 49 ' 


16.90" N 
52.70" W 


SCI = +0.36233281 
SC 2 = -0.36296869 
SC3 = -0.35162108 


D 


N. Dakota 


46° 21* 
98° 20' 


57.20" N 
08.77" W 


SD 1 = +0.72144222 
SD2 = -0.1 0039077 
SD3 = -0.68515898 


E 


LaReunion 


20° 58* 


26.47" S 


SE1 = -0.35585299 




Island 


55° 17' 


24.25" E 


SE2 = +0.53214375 
SE3 = +0.76823587 


P 


Argentina 


43° 03' 
65° 11 ' 


12.53" S 
27.60" W 


SF1 = -0.68022941 
SF2 = +0.30756230 
SF3 = -0.66535207 


G 


Trinidad 


10° 42' 
61° 38* 


06.20" N 
20.30" W 


SGI = +0 . 18449565 
SG2 = +0.46687116 
SG3 = -0.86486570 


H 


Japan 


340 451 
129° 27 • 


53.26" N 
12.49" E 


SHI = +0.56547261 
SH2 = -0.52409936 
SH3 = +0.63633639 



TA3LE I. OMEGA Stations 
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FIGURE I - OMEGA SYSTEM TRANSMITTER LOCATIONS 
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FIGURE 2 - OMEGA SYSTEM TRANSMISSION PATTERN 
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transmitted signal appears as a standard sinusoidal wave. 
Assuming no attenuation of the signal, a perfect waveguide 
formed by the earth and the ionosphere for the signal to 
propagate through, no interfering carrier or noise signal, 
and the transmission traveling at the speed of light, the 
received signal would appear as a sinusoidal wave but offset 
at the receiver in time t by d/c, the time necessary for the 
signal to travel distance d at the velocity of light c. 




However, this ideal situation is rarely achieved. The 
propagated signal is affected by several factors described 
in a subsequent section. Due to the repetitive nature of 
sinusoidal waves, the precise number of wave lengths must be 
computed between the transmitter and the receiver. Also, 
the absolute time t of the transmitter must be known. These 
variables are calculated in the combinational filter. 



Position fixing with the use of OMEGA signals is 
accomplished by two methods. The first method utilizes 
phase difference information between stations and the 
resultant nypermolic lines of position. The second method 
uses circular lines of position obtained by measuring phase 
with respect to a stable frequency standard. The P-3C Omega 
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navigation system uses the circular lines of position or 
rho-rho method to determine position fixes. 

Of primary importance in utilizing the OMEGA VL? signals 
for the rho-rho method are the ability to detect and measure 
the phase of a received signal, and the ability to predict 
the phase of the signal. In order to detect and measure the 
phase of a received signal using the rho-rho method, the 
first process performed in the Omega navigation system is to 
relate the internal clock, of the computer with the 
transmission burst pattern of the stations. This 
synchronization process is described in more detail in a 
subsequent section. The overall effect allows the computer 
program to know what station and frequency is being 
processed at any internal system time t. If all three 
frequencies from all eight stations could be received at any 
point on the earth's surface, the ten second reception 
pattern would appear as in Figure 3. For a more realistic 
reception pattern, Figure 4 depicts the three closest 
transmitter stations to Hawaii, and Figure 5 represents the 
signal pattern an aircraft would receive at Hawaii. 

Once synchronization is completed, the Omega navigation 
system commences to process the received signals through 
various software filters in order to measure the phase. The 
filters which process the received signal are composed of 
three burst filters (ie. one for each frequency) and 24 
tracking filters (ie. one for each of the three frequencies 
on each of the eight stations) . The primary function of the 
burst filters is to calculate the phase measurement and the 
phase variance for each frequency and station received by 
the receiver. The tracking filter, on receiving the phase 
measurement and variance of the phase measurement from the 
burst filter, utilizes these values to acquire a more 
accurate estimate of phase and phase variance over several 
t en second time intervals. 
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FIGURE 3 -OMEGA SYSTEM RECEIVER PATTERN 
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FIGURE 4 - SPECIFIC TRANSMITTER LOCATIONS 
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FIGURE 5 - SPECIFIC RECEPTION PATTERN 



The tracking filter performs two different updates on 
the burst filter data. The first update is a time update. 
This update is required because of the displacement of the 
aircraft receiver between successive measurements of the 
received signal. The time update procedure projects the 
last estimate of phase through time so that it reflects the 
new position of the aircraft. These procedures are defined 
as rate aiding and are dependent on velocity sources 
external to the Omega system. In estimating the change in 
position, it is possible to estimate what the change should 
be for a particular station and frequency. This predicted 
phase change is then added to the last estimate of phase to 
create a new phase estimate. 
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measurement update, and is perfor 
filter supplies a measurement of 
the tracking filter. This mea 
combines the rate aiding phase 
filter with the burst filter phas 
better estimate of the phase 
equation . 



(t) + A 'PEER (t-1) ] ’At 

rate which represents the 
along the arc connecting a 
pERR(t-l) is the estimate of 
and is calculated in the 
The second update is a 
med every time the burst 
phase and phase variance to 
surement update procedure 
estimate from the tracking 
e measurement to produce a 
based on the following 



p = pJK - PTRK(t) + 



A0DR(t) + A PEER ( t ) ] * (0.425) 



where p is defined to be the angular difference of phase. 
The combination of the phase measurement from the burst 
filter and the phase estimate from the tracking filter is a 
weighted average, where the coefficients are computed using 
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the phase variance from the burst filter (o~ 2 (pJK) and the 
phase variance from the tracking filter (cr z 0TRK) . 



07RK (t) 



<DTRK (t-1) 



+ ARCTAN 



(cr 2<f)TRK * S) 

( (cr 2 0JK + cr 2 0T RK) 




where S = SIN <J» and C = COS 0. 



Once the measurement of phase has been calculated, it is 
supplied to the combinational filter for processing along 
with a phase estimate calibrated in the propagation 
prediction procedure. With the measured phase and the 
predicted phase available for position fixing, the last 
requirement to be implemented for utilization of the rho-rho 
technique of position fixing is to synchronize the receiver 
oscillator with the transmitter oscillator. This is 
accomplished in the combinational filter by computing the 
time difference between the transmitting station phase and 
the receiver oscillator phase. 



1 • Signal Coverage 

The VLF signals utilized by the OMEGA system are 
transmitted ever extremely long distances of the earth's 
surface, because they are propagated through a natural wave 
guide formed by the surface of the earth and the ionosphere. 
Due to tne advantageous properties of high phase stability 
and low attenuation rates of VLF signals, no modulation of 
the signals is required prior to transmission. Because these 
signals are unmodulated, only their relative position in any 
ten second time interval enables receivers to identify the 
station being received. Df primary importance in utilizing 
the OMEGA VLF signals are the ability to detect and measure 
the phase of a received signal, and the ability to predict 
the phase of the signal. Several factors have an effect on 



the propagation of VLF signals and consequently the ability 
of an airborne Omega system to receive and process the 
signals fcr accurate positioning [Ref. 6]. These factors 
are described in the section pertaining to propagation 
prediction. 



2 • Rece iver -Converter 



The AN/ARN-99 receiver was developed by the 
Electronics Division of the Northrop Corporation in 
Havthrons, California. It was designed for operation with 
all eight transmitting stations and to provide a continuous 
update of aircraft position. The receiver- con verter 
installed in the P-3C consists of several modular assemblies 
as shown in Figure 6. The antenna, composed of two loop 
antennas fixed at right angles to each other and mounted at 
45° angles to the aircraft center line, is controlled by the 
burst filter function in the computer through the antenna 
switching matrix. The three heterodyne receivers have RF 
and IF sections for processing of the input signals. The 
requirements for filtering the received signals in these 
sections are very stringent. The filtering in the RF 
section removes the RF image, including all frequencies in 
the area of the harmonics of the local oscillator. The 
narrow band filtering in the IF section provides for the 
rejection of interfering carriers. In addition, limiters 
are used for controlling the dynamic signal levels, and a 
bandpass filter is used to remove the harmonics created by 
the previous limiters. The resultant output of the RF and 
IF filtering sections leaves the fundamental frequency as a 
sinusoid. This output is fed into a sine and cosine 
correlator which utilizes a square reference signal to 
produce an output of two direct current signals. One dc 
signal is proportional to the sine of the burst phase, and 
the other dc signal is proportional to the cosine of the 
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FIGURE 6 - RECEIVER-CONVERTER/COMPUTER INTERFACE 



III. 



IN IEEN AL PR OC ESSING FUNC TIONS 



The Omega navigation systeQ in the P-3C aircraft is one 
of frve navigation subsystems which comprise the avionics 
navigation system. The remaining subsystems, in their 
priority of use, are a primary inertial, a secondary 
inertial, a doppler radar, and an air data computer. Tne 
selection and utilization of each subsystem is dependent on 
the control of tne Navigator/Communicator and on the 
availability of the subsystem. A block diagram of the 
navigation subsystems is shown in Figure 7. 

A. OMEGA PROCESSING 

The internal processing functions of the Omega 
navigation system are shown in Figure 8 [Ref. 4], A more 
detailed description of the internal processing functions is 
presented in Figure 9 showing the pertinent data transferred 
between the functions. 



1 . Initialization 



Initialization of the Omega navigation system begins 
with the insertion of the date and time entries by the 
Navigator/Communicator and subsequent activation of the 
initial cn top condition or mark aircraft position. The 
date and time entries are used within the propagation 
preai ction module to calculate an estimate of the phase 
value derived from current aircraft position at the moment 
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FIGURE 7 - NAVIGATION SYSTEM BLOCK DIAGRAM 




FIGURE 8 - INTERNAL OMEGA PROCESSING FUNCTIONS 
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FIGURE 9- DATA FLOW GRAPH OF INTERNAL FUNCTIONS 





of initialization. These phase values are used to 
initialize the tracking filters. The variance values within 
the combinational filter are also initialized at this time, 
and the position matrix R, utilized in the velocity and 
navigation processing function, is initialized based on 
current aircraft latitude and longitude. 

2 • Hardw are Testing 



Hardware testing begins immediately after the 
termination of the initialization function, and also after 
each restart. The main objective of the hardware testing 
module is to determine the status of the Omega 
receiver/converter. Additionally, hardware testing provides 
a method of assuring proper operation of all Omega 
equipment. The testing involves a series of subtests which 
validates the operation of the input/output interface 
between the receiver/converter and the computer. The 
subtests, which are described in greater detail in Ref. 4, 
are a communication subtest, a coherence status subtest, an 
Omega output subtest, an Omega input subtest, a phase angle 
to digital subtest, a phase counter subtest, and a RF/IF 
subtest. An oscillator drift subtest, which is classified 
as a hardware test, is not implemented until single 
frequency processing commences in the combinational filter 
routine . 



3* Synchroniz at ion 

Synchronization of the Omega navigation system to 
me Omega transmission pattern is the most critical function 
executed by the navigation system. 3y aligning itself with 
the ten second interval transmission pattern, the Omega 
navigation set processes the correct VLF signal bursts in 
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subsequent routines. Synchronization is accomplished by 
utilizing the variable length transmission time intervals of 
the various frequencies and a differential correlation 
technique. The first assumption made in the routine is 
that all eight stations that are presently in operation are 
operating on all three frequencies. The routine, starting 
at some random time, accepts from the receiver-converter a 
tec second interval of data composed of 100 sine and 100 
cosine values over contiguous 100 millisecond intervals. 
One hundred correlation coefficients are calculated for the 
frequency under consideration utilizing the difference 
between short sum and long sum intervals and adding these 
differences over the eight count burst pattern. Figure 10 
depicts the short sum and long sum intervals over a ten 
second interval pattern, and the limits of the short and 
long sums. The correlation coefficient formula is given 
below . 



8 


,2 


/ I ^ U \ 


2 / I+JU+2 >2 


1 \ 


2 


c <d = 2 


2 x (s) * 1 


^ Y(K)' 


|- X x(K > H 


X r ‘ s > 




j=i 


\K= 1+ JL / 


U=I+JL j 


' \K=I + JL-2 / 


U=I + JL-2 j 




(1=1, 100) 


JU-JL+1 


JU-JL+5 






The values X (K) and Y (K) are the sine and cosine values 
respectively received from the receiver-converter. Each of 
the confidence coefficients can be considered to be a 
measurement cf how well the received ten second burst 
pattern matches the transmitted burst pattern. The higher 
the coefficient value, the greater confidence is associated 
to the index number I being the starting point of the 
transmitted burst pattern within ± 50 milliseconds. The 

four largest confidence coefficients are saved in descending 
order of magnitude, along with their respective index number 
I. This index number is used in determining the starting 
point of the ten second burst pattern. Various conditions 
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relating to the four highest confidence coefficients and 
index nuafcers are tested in order to determine the success 
or the failure of the synchronization attempt. If the 
current confidence level fails, the next sequential 
frequency is utilized to acquire the next ten second 
interval cr data. The synchronization routine continues in 
this mode until a successful synchronization occurs or until 
the Navigat or/Comm unicator intervenes. 
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FIGURE 10- SYNCHRONIZATION PARAMETERS 
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are computed to improve the accuracy of the phase 
measurement., and for the issuing of antenna selection and 
test selection commands. Reference 4, pages 86 to 101, 
describes in detail the computations involved in calculating 
the phase and the phase variance for each station and 
frequency used by the Omega navigation system. 
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5 • Trac king Filters 



Once the burst filter processing has acquired a 
phase measurement and an estimate of the variance of the 
phase measurement for a specific station and frequency, 
these values are used as inputs to only one of the 24 
tracking filters. Each filter is similar in the processing 
of the data, tut each is different in its computational 
constants for the various frequencies and stations. A 
simplified description of a tracking filter is that it is 
utilized tc acquire a better estimate of the phase 
measurement by averaging the outputs of the burst filter 
over successive ten second intervals. The tracking filter 
routine, using the measured phase from the burst filter and 
an updated estimate of phase based upon previous 
measurements and dead reckoning information, combines the 
measured and estimated phase and utilizes a weighted average 
to compute a new phase estimate and phase variance estimate. 
Also computed within the tracking filter is the variance of 
the estimate of the error in the phase rate of change due to 
aircraft displacement in time. Once these three values are 
computed for a particular station and frequency, a test is 
made to determine the accuracy of the phase measurement. If 
the variance of the phase measurement is less than or equal 
to (0.06 it radians) 2 , the phase measurement is considered 
accurate enough for computations within the combinational 
filter. After three successive measurements and successful 
tests, a flag is set indicating that the data is valid and 
available, and the combinational filter reads the outputs 
from the tracking filter. Once the outputs have been read, 
the count of successive measurements is reset to zero, the 
variances of the phase measurements are set to the 
initialization values, and the data valid and available flag 
is cleared. If at any time in the tracking filter the 
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variance cf the phase measurement exceeds (0.06 tt radians) 2 , 
the count of the successive measurements is reset to zero, 
and the data valid and available flag is cleared. 



6 . Combinational Fil ter 

The input values to the combinational filter from 
the tracking filter are the refined estimates of phase 
measurement, 0TRK, phase variance, cr 2 <J)TRK, and phase rate 
variance, cr 2 <l)E2R. The inputs from the propagation 

prediction routine are the estimates of phase, <j!PR0P, 
improved for propagation effects and phase variance, 
a" 2 <J>PR0P. Within the combinational filter routine, these 
input values are statistically combined to provide the best 
estimates of system position and velocity. The 
combinational filter routine provides for the conversion 
from phase tc geodetic coordinates (latitude and longitude) , 
and also determines the phase and frequency differences 
between the receiver oscillator and the transmitted OMEGA 
signal. Cnee the Omega receiver oscillator is calibrated, 
the Omega navigation system operates when only two 
transmitting stations are accessible. In addition, the 
combinational filter is used for lane determination through 
the technique of multiple state vectors. 



a. Description 



The combinational filter is a Kalman filter. 
Reference 4 describes the Kalman filter technique to 
filtering and prediction as a linear, recursive, minimum 
variance filter. R. G. Brown and L. L. Hagerman [Ref. 9] 
describe a Kalman filter as simply a means of estimating the 
various states of a random process from a set of discrete 
measurements having a known linear connection to these 
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states. It was further noted in Ref. 10, that a Kalman 
filter operates only on the system errors and not on total 
quantities such as position and velocity. The basic concepts 
involved are those of stare, state transition, measurement, 
and optimal weighting [Ref. 11]. The states of the 
combinational filter are differentials of system variables 
[ie. error in position (<5<J)2 and 6<b3) , oscillator start time 
(t 0 ) , rate of change in oscillator start time ( t 0 ) , and 
error in velocity (&VC2 and <$VC3) ]. The following chart 
describes the elements of the state vector X [6 x 1]. 



ELEMENT SYMBOLS 

1 6<t 2 

2 6(1)3 

3 t 0 



4 t c 

5 Sv C2 



6 SVC 3 



MEANING 

Position error along the -R2 direction of 
the system position matrix R. 

Position error along the -R3 direction of 
the system position matrix R. 

Time difference between the transmitting 
station phase and the receiver oscillator 
phase. 

Time rate of change of t c . 

Error in the east component of velocity 
which is error along the R2 direction of 
the system position matrix R. 

Error in the north component of velocity 
which is error along the R3 direction of 
the system position matrix R. 



The state of the system is described by the 
solution of linear vector differential equations depicting 
system error growth. The Kalman filter method linearly 
combines the previous estimate of the state vector, X, with 
a measurement, to approach a minimum variance estimation. 
This minimum variance estimate is then time updated until 
the following measurement. The measurement calculation is 
explained by the following equations. 
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Y = <J>TRK - 0PROP 



where Y is the difference between the tracking filter phase 
measurement and the propagation prediction phase estimate of 
the phase measurement. 

Residual = Y - M*X(t-1) 

where Residual (RES) is the difference between the 
calculated error in the measurement, Y, and the system 
estimate of the error, H*X (t-1) . The state vector X(t-1) is 
the previous estimate of the state vector, and a [ 1 x 6] is 
the measurement matrix described below. 



ELEMENT 



1 



2 



3 



EQUATIONS 

3 



EMER * 




R2J*SIJ 



J= 1 

K * SIN(DD) 



EMSR 



*1 

i— i 



K * 



R 3 J*SI J 



SIN (<DD) 



VL 

AK 



MEANING 

EMER is the earths mean equa- 
torial radius = 20925741.47 ft 
R2J and R3J are the elements of 
the R2 and R3 vectors of the 
position matrix R. 

SIJ’s (I=A,3,..,H) are the 
elements of the station 
direction cosine values. 

VL is the average propagation 
velocity, and AK (K=1,2,3) is 
the wave lengths of the 
frequencies 10.28,11.3, 13.6kHz 



where <t)D 



3 



ARCCOS 



<2 



R 1J*SIJ) 



and the remaining elements of 



M are zero. New, minimum variance, state vectors are 
produced by using the following equations. 
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Xnew = Xold + b * (RES) 



where b [6 x 1 ] is the optimal weighting vector matrix. 

This optimal weighting vector, b, is established 
by means of time propagation of system error growth in 
combination with information contained within the 
measurement. The weighting is a function of the covariance 
matrix, P [6 x 6], of the difference between system error 
state vector, X, and the average state vector Xavg. 

P = E * [ (X - Xavg) * (X - Xavg) T ] 

where E is the expected value of the difference between the 
value of the difference between the value for the real stare 
vector X, if all conditions were known, and the average 
state vector, Xavg. The following chart describes the 
elements of the covariance matrix P. 



ELEMENT 


SYMEOLS 


MEANING 


PI 1 


cr 


2<t>2 


The position variance for the position 
error, <$<J>2, along the -R2 direction 


P22 


cr 


2<D3 


The position variance for the position 
error, <$<i>3 , along the -R3 direction 


P33 


c r 


O 

(\l 


The variance of oscillator phase offset 


P44 


cr 


2 t c 


The variance of oscillator drift rate 


?55 


cr 


2VC2 


The variance of the error in the east 
component of velocity 6v C2 


P66 


cr 


2 VC3 


The variance of the error in the north 
component of velocity <SvC3 



All remaining elemenrs of the covariance matrix are zero. 

The combinational filter effectively computes 
optimal estimates cf position and velocity along with the 
remaining elements of a state vector. The filter uses an 



37 



observation of the difference between the phase of the 
tracking filter and a value of the phase formulated on the 
time updated position of the combinational filter. Based on 
this measurement, the combinational filter derives estimates 
of position errors, velocity errors, oscillator start time 
and oscillator drift. 

Eecause the Omega navigation system functions 
with several distinct velocity sensors (ie. inertial, 
doppler, or air data) , a different mode of operation exists 
for the combinational filter for each- sensor. Each method 
requires the filter to estimate, predict, and control a 
different set of system errors. Each set of system errors 
conforms to the dead reckoning velocity source used for time 
updating . 



The combination of uncertainties in position 
error and receiver oscillator start time produces a 
vagueness of the number of wave lengths or lanes between the 
transmitting station and the receiver. This ambiguity occurs 
primarily when the system is first initialized. To find the 
solution to this lane ambiguity, several estimates of the 
stare vector, X(i) 's, are calculated by the filter. These 
X(i)'s correspond to the different possible lanes or 
integral values of phase measurement. As measurements are 
acquired, corresponding to the three frequencies, the 
uncertainity in oscillator start time, to, is diminished. 
Measurements from different stations tend to minimize the 
position uncertainity. The general effect is a reduction in 
the number of state vector estimates, X (i) 's, that can be 
considered logical estimates of the most accurate state 
vector until only one credible estimate remains. The 
criterion for reasonableness is based on variance 
considerations. 
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b. Operation 



There are three distinct operations identified 
with the combinational filter routine. 

1) Initialization 

2) Time update 

3) Measurement update 



The initialization operation assigns the initial 
values of variance to the covariance matrix P, sets the 
elements of the initial state vector to zero, and sets the 
system driving noises as a function of the velocity source. 

In the time update operation, both the state 
vector X, and the covariance matrix P, are predicted from 
the last update. The following equations are utilized in 
the time update operation. 

Xnew = |new * Xold 

Pnew = |new * Pold * (jjnew T + Nnew 



The transition matrix <j[ [ 6 x 6] mathematically expresses the 
propagation of errors across the time interval A t since the 
last update. The elements of the transition matrix are 
shewn in the following chart. 



1 

0 

0 

0 

0 

0 



0 

1 

0 

0 

0 

0 



0 

0 

1 

0 

0 

0 



0 

0 

At 

1 

0 

0 



At 

0 

0 

0 

1 - £i*At 
0 



0 

At 

0 

0 

0 

1-pi*At 



where pi (i=I,D, or A) signifies the inverse correlation 
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time of &VC2 and &VC3 for the velocity source mode utilized. 
The transition matrix $) furnishes the propagation of the 
predicted elements of the state vector X and the covariance 
matrix P across the time interval A t. The variance effects 
across this time interval are propagated by the additive 
diagonal noise matrix N [6 x 6], The following chart depicts 
the elements of the noise matrix N. 



0 

0 

0 

0 

0 

0 



0 

0 

0 

0 

0 

0 



0 

0 

0 

0 

0 

0 



0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 (5i*cr 2 i*At 0 

0 0 £i*o" 2 i*At 



where (J i is the inverse correlation time described 
previously, and c r 2 i is the variance of the velocity errors 
(<S V C 2 and <5 V C 3) and is dependent on velocity source utilized 



(i=I , D , or A) . 



During the time update operation, the multiple 
state vectors, X (i) 1 s, which exist until the correct lane 
has been determined must also be time updated. The time 
update equation for the state vector X must be sequenced 
through all state vectors. 

X (i) new = |new * X (i) old for all i 

The covariance matrix time update equation is computed once, 
because there is only one covariance matrix regardless of 
the number of state vectors. 



After completion of the time update operation 
and prior tc the measurement update operation, the 
measurement residual, (RES) , measurement matrix, M , and the 
measurement confidence scalar, C, are calculated for each 
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tracking filter input. The measurement residual and 
measurement matrix calculations were previously explained in 
the section describing the combinational filter. The 
measurement confidence scalar, C, is a measurement of noise 
computed by the following equation. 



C = <r 2 $prop + cr 2(J>TRK 



where cr 2 (£prop and cr 20TRK are the variance of 
the propagation prediction phase estimate and the tracking 
filter phase measurement. 



The measurement update operation 
generation of a linear, unbiased weighting 
computation cf a new, minimum variance, state 
the formula listed below. 



involves the 
vector b, and 
vector X by 



Xnew = Xold + b * (RES) 

To generate the optimum weighting vector, b, three 
expressions are reguired. 

1) the predicted measurement variance - H*P*M T 

2) the measurement confidence scalar - C 

3) the divergence control factor - € (epsilon) 

The weighting vector is calculated from the following 
equations. 

Q = M*P*M T + C 
b = (P*M T + €TM T ) / Q 



where Q is defined to be a scalar quantity depicting the 
lane variance, and e is a factor to prevent divergence of 
the estimation process. Also updated in the measurement 
update process is the covariance matrix based on the 
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following formula. 



Pnew = Pold - b * M * Pold + e * m T * b T 

The retention of new state vectors is predicated on the 
testing cf the position elements in the state vector, and 
the variance values in the covariance matrix. 

7 • Pr op a gati on Predic t ion 

In order to effectively resolve the aircraft 

position from the phase information received, the phase 

velocity cf the wave along the propagated path and the 

length of the path must be known. Several factors or 

♦ 

effects create imperfections in the wave guide through which 
the VLF signals propagate. The propagation prediction 
function establishes the best estimate of phase and phase 
variance taking into account the factors listed below which 
effect the phase velocity of the propagated signals. 

a. Diurnal Effects 

The positional changes of the sun over the earth 
adjust the size and the shape of the ionosphere. The change 
of the wave guide through which the VLF signals propagate 
causes a change in the velocities of the phase signal which 
results in a difference in phase angle received at the Omega 
receiver location. 

b. Ground Conductivity 

This factor takes into account the effect of 
various mediums over which VLF signals travel. Water, which 
is a near perfect conductor of VLF signals, does not greatly 
affect the signals. Land, which is a less perfect conductor 
of VLF signals, does affect the signals to a greater extent. 
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Conductivity patterns have been measured and are predictable 
[Ref. 6]. 



c. Earth's Magnetic Field 

This factor affects both the attenuation rate 
and the velocity of the signal. The phase shift on 
refraction of the VLF signal off of the ionosphere depends 
on the interaction of the signal with the electrons, and 
this depends on the orientation of the magnetic field with 
respect to the direction of propagation and on the magnitude 
of the field [Ref. 12]. 

d. Latitude and Spheroidal Effects 

This factor takes into account the nonspherical 
shape of the earth and the adjusted path of travel from the 
transmitting station to the Omega receiver. 

e. Folar Cap Absorption (PCA) 

VLF signals which pass over the Polar Cap 
experience an abnormal rate of absorption during both day 
and night. PCA produces large changes in VLF signal patterns 
which can last for several days. 

f. Sudden Ionospheric Disturbance (SID) 

Sclar flares emanating on the sun's surface 
increase the ionization within the atmosphere. These SID's 
cause large changes in the phase of VLF signals passing 
through the area of activity. 

Of the six factors which affect the propagation of 
VLF signals, the first four are predictable and can be 
compensated for algorithmically, while the last twc are 
unpredictable and cannot be compensated for within an Omega 
navigation system. 
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IV. VELOCITY AND NAV IGAT ION PROCESSING FUNCTION 



The velocity and navigation processing function is 
divided into two routines in the Omega navigation system. 
The velocity processing routine, utilizing the current 
navigational source (ie. inertial, doppler, or air data), 
calculat-es the velocities along the system axes from the 
available source velocities. These system axes velocities, 
combined with the state vector corrections from the 
combinational filter, are supplied to the navigation routine 
in order to update the Omega system position matrix and the 
Omega aircraft latitude and longitude position. 

A. VELOCITY PROCESSING 

In the velocity processing routine, the aircraft current 
true airspeed velocity (VTAS) and heading (ACHDG) , the 
current wind direction (WD) and velocity (VW) , and the 
dopplers current velocity components along heading (VDA) and 
across heading (VDC) are utilized to compute the velocities 
in the North/South and East/West directions from the 
following equations, 

VDN = VDA * CDS (ACHDG) - VDC * SIN (ACHDG) 

VDE = VDC * COS (ACHDG) + VDA * SIN (ACHDG) 

VAN = VTAS* COS (ACHDG) + VW * COS (WD) 

VAE = VTAS* SIN (ACHDG) + VW * SIN(WD) 

where VDN and VDE signify doppler velocities North/South and 
East/West, and VAN and VAE signify air data velocities 
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North/South and East/West respectively. In constructing the 
data flow graph of these equations, it was self-evident that 
some variables could be calculated, stored in memory, and 
retrieved from memory when they were required in later 
calculations. Figure 12 depicts the data flow graph for the 
North/South and East/West computations. Implementing PL/M 
code directly from this data flow graph was straight 
forward. 

On completion of the North/South and East/West 
velocities computations, a test is made in the velocity 
processing function to determine if the navigation mode has 
changed since the last iteration of the velocity function. 
If an upgrade of the navigation mode has occurred (ie. 
doppler to inertial, or air data to doppler or inertial) , 
then the summations (DELVC2 and DELVC3) of the corrections 
to the system velocity components from the combinational 
filter are set equal to zero, and new velocities along the 
system axes are calculated. Figure 13 depicts the 
navigation mode test and applicable equations for computing 
the system axes velocities. 

In order to save execution time in the calculation of 
the formulas and to simplify the PL/M code, the functions 

(ViE * COS (SHDG A) + ViN * SIN {SHDGA) ) 

(ViN * COS (SHDGA) - ViE * SIN (SHDGA) ) 

were calculated prior to the navigation mode test. Figure 14 
depicts the data flow graph for these two functions. After 
successful computation of the above functions and the 
navigation mode test, the corrected velocities (VC2 and VC3) 
along the system axes are calculated and the resultant 
values used in the navigation routine. 
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FIGURE 12 - NORTH/SOUTH/EAST/WEST VELOCITY DATA FLOW GRAPH 
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FIGURE 13 - VELOCITY PROCESSING FLOWCHART 



COS(SHDGA) 
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FIGURE 14 - PARTIAL FUNCTION DATA FLOW GRAPH 



3 . 



NAVIGATION PROCESSING 



In the navigation processing routine, the Omega system 
position matrix R [3 X 3] is initialized, the angular 
rotations (DELT2 and DSLT3) about the . system axes are 
calculated using the corrected velocities supplied by the 
velocity processing routine and the state vector corrections 
from the combinational filter, the system position matrix is 
updated, and a new Omega aircraft latitude (OLATAC) and 
longitude (GLCNAC) are calculated. The flowchart of the 
navigation routine is shown in Figure 15. 

The initialization of the Omega system position matrix R 
occurs during the initialization function when the 
Navigator/Communicator inserts the aircraft initial latitude 
(OLATIN) and longitude (OLONIN) . These values, with an 
initial condition that the system heading angle (SHDGA) is 



zero , 


are input to 


the equations listed below 


R 1 1 = 


SIN (CIATIN) 






R 1 2 = 


CCS (OLATIN) 


* 


COS (OLONIN) 


R 1 3 = 


COS (OLATIN) 


* 


SIN (OLONIN) 


E21 = 


COS (OLATIN) 


* 


SIN (SHDGA) = 0 


R22 = 


-SIN (OLONIN) 


* 


COS (SHDGA) - 




SIN (OLATIN) 


£ 


COS (OLONIN) * SIN (SHDGA) 


= 


-SIN (CLCNIN) 






R23 = 


COS (OLONIN) 


a, 

*v 


COS (SHDGA) - 




SIN (OLATIN) 


* 


SIN (OLONIN) * SIN (SHDGA) 


= 


COS (CLCNIN) 






R31 = 


COS (OLATIN) 


* 


COS (SHDGA) = COS (OLATIN) 
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R32 = SIN (OLONIN) * SIN(SHDGA) - 

SIN {OLATIN) * COS (OLONIN) * COS(SHDGA) 
= -SIN (CIATIN) * COS (OLONIN) 

333 = -COS (OLONIN) * SIN ( SHDGA) - 

SIN (CIATIN) * SIN (OLONIN) * COS (SHDGA) 
= -SIN (OLATIN) * SIN(OLONIN) 



These nine eguations represent nine direction cosine values. 
These values represent the orientation of each of the three 
system axes (31, 32, and 33) in the earth's fixed coordinate 
frame. In constructing the data flow graph for the 
initialization equations, the concept of minimizing 
execution time prevailed, and as variables were calculated, 
they were stored in memory and retrieved as they were 
required in later calculations. Figure 16 portrays the data 
flew graph for the initialization of the position matrix 3. 
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FIGURE 15 - NAVIGATION PROCESSING FLOWCHART 
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R I I 




FIGURE 16 - INITIALIZATION OF POSITION MATRIX R 
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